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To Percy Deift at the occasion of his 60 'th birthday 
Abstract 

In this paper we compute the Birkhoff normal form of the periodic 
Toda lattice up to order four. As an application, we verify that Kol- 
mogorov's nondegeneracy condition in the KAM theorem holds almost 
everywhere in phase space. 1 

1 Introduction 

Consider the periodic Toda lattice with period N (N > 2), 

q n = d Pn H, p n = -d qn H 

for n G Z, where the (real) coordinates {q n ,Pn)nez satisfy (q n +N ,Ph+n) = 
(<ZnjPn) for any n € Z and the Hamiltonian H is given by 
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with potential V{x) = j 2 e Sx + V x x + V 2 and 7, 8, Vi, V 2 £ K constants. The 
Toda lattice has been introduced by Toda ^2] an d studied extensively in the 
sequel. It is an FPU lattice, i.e. a Hamiltonian system of particles in one space 
dimension with nearest neighbor interaction. Models of this type have been 
studied by Fermi-Pasta-Ulam [FPU]. In numerical experiments they found re- 
current features for the lattices they considered. Despite an enormous effort 
from the physics and mathematics community, some of these numerical experi- 
ments still defy an explanation. For a recent account of the fascinating history 
of the FPU problem, see e.g. pQ. At least in the case of the periodic Toda 

* Supported in part by the Swiss National Science Foundation 

^Supported in part by the Swiss National Science Foundation, the programme SPECT and 
the European Community through the FP6 Marie Curie RTN ENIGMA (MRTN-CT-2004- 
5652) 

1 2000 Mathematics Subject Classification: 37J35, 37J40, 70H08 



1 



2 



1 INTRODUCTION 



lattice, the recurrent features can be fully accounted for. In fact, Flaschka [2], 
Henon [3], and Manakov |S] independently proved that the periodic Toda lattice 
is integrable. In this paper, we show that (small) Hamiltonian perturbations of 
the Toda lattice have many - large and small - quasi-periodic solutions as well. 
To this end, we put the periodic Toda lattice into Birkhoff normal form up to 
order 4 and then show that Kolmogorov's nondegeneracy condition of the KAM 
theorem holds almost everywhere in phase space - see e.g. 0, Appendix G, for 
the notion of Birkhoff normal form up to order m. 

Returning to the Toda Hamiltonian H, let us first note that w.l.o.g. we 
can assume that V% = V2 = 0. When expressed in the canonical coordinates 
(Sqj, jPj)i<j<N the Hamiltonian H is, up to a scaling factor S~ 2 , of the form 

N N 

H Tod a=2Y,Pn+ a2 Y, eqn ~ qn+1 > C 1 ) 
n— 1 n—1 

where a 2 — (jS) 2 . Moreover notice that the total momentum J2n=iP n * s con " 
served; hence the motion of the center of mass jj X) n =i Qn 1S linear and therefore 
unbounded. However, it turns out that the orbits of the system relative to the 
center of mass all lie on tori, generically of dimension N — 1. To describe 
these orbits, it suffices to consider the relative coordinates v n := q n +i — Qn 
(1 < n < N — 1) and their conjugate ones, u n := n/3 — Y%=i Vk (1 < n < N— 1), 
where f3 = jj ^2f=iPn- The corresponding phase space is then R 2Ar ~ 2 . 
In terms of the variables (vk, Uk)i<k<N-i, Hroda is given by 

H,, a = ^f+l U/^ut-u^f+u^+a 2 (j^e-^+e^A , (2) 

with parameters a and f3 - see section [21 for more details. 
The main result of this paper is the following 

Theorem 1.1. For any fixed (3 G K, a > 0, and N > 2, the periodic Toda lattice 
admits a Birkhoff normal form. More precisely, there are (globally defined) 
canonical coordinates (xk,yk)i<k<N-i so that Hp_ a , when expressed in these 
coordinates, takes the form l Kp. a (I) := -^y — h H a (I), where H a (I) is a real 
analytic function of the action variables Ik = (x 2 + yl)/2 (I < k < N — 1). 
Moreover, near I = 0, H a (I) has an expansion of the form 

N-i ^ N-l 

Na 2 +aY; Skh + E 7 ' + °( /3 )' ( 3 ) 

fe=l k=l 

with s k = 2sin ^ for 1 < k < N - 1. 

Remark 1.2. In particular, Theorem savs that near 0, the Toda lattice 
can be viewed as a system of (N — 1) nonlinear harmonic oscillators which are 
uncoupled up to order 4. The system is resonant at with resonance lattice 
Res := {(fcj^KjXjv-i £ Z N ^ X \ X)j=i ^j s j = 0} °f dimension at least i_ jV 2 ~ 1 j as 
Sk — SN-k = for any 1 < k < l^^j. 
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Corollary 1.3. Let a > and (3 eK be arbitrary. Then the Hessian ofKp^ a {I) 
at I = is given by 

In particular, the frequency map I i— > V/5{^ iCt is nondegenerate at I — and 
hence, by analyticity, nondegenerate on an open dense subset of (]R>o) Ar_1 - 

Theorem 11.11 and Corollary 11.31 allow to apply the KAM-theorem (cf e.g. 
[7| to the Toda lattice on an open dense subset of the phase space. Moreover, 
as the Hessian of JCp, a at I = is positive definite and hence Jip,a is strictly 
convex near 1 = 0, one can also apply Nekhoroshev's theorem (cf e.g. ^0]) to 
the Toda lattice near the equilibirum point. 

Using the method of proof of this paper we computed in [S] the Birkhoff 
normal form up to order 4 for any FPU chain and applied these computations 
to improve on results of Rink 11 with regard to the verification of Kolmogorov's 
condition near the equilibrium point for these systems. 

Acknowledgement: It is a great pleasure to thank Percy Deift and Yves Colin 
de Verdiere for helpful discussions. 



2 Relative coordinates 

To prove the integrability of the Toda lattice, Flaschka introduced in |2] the 
(noncanonical) coordinates 

b n :=- Pn eR, a n := ae i(?»-««+0 g R >0 ( n € Z). (4) 

They are coordinates which describe the motion of the Toda lattice relative 
to the center of mass. In these coordinates, the Hamiltonian Hroda takes the 
simple form H = i Xm=i + Sn=i a «' ana - ^he ecma tions of motion are 

Note that (b n +N, clu+n) = {b n , a>n) for any n € Z, and O^Li °« = aN ■ Hence we 
can identify the sequences (6 n )nez and (a n ) n ez with the vectors (6„)i<, i <Ar 6 
K w and {a n )i< n <N € K>o- The phase space of the system © is then given by 

M := R N x R* , 

and it turns out that (jSJ) is a Hamiltonian system with respect to the non- 
standard Poisson structure J = J{, a , defined at a point (b, a) = {b n ,a n )i< n <N 

by 

A 
-A T 



J 
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where A is the 6-independent N x TV-matrix 

(a x ... 



A = l 
2 



-Oil 0-2 
—da «3 



-CLN \ 







a N J 



(6) 



\ ... -a^! 

The Poisson bracket corresponding to J is then given by 

{F,G}, } {b,a) = ((V b F,V a F),J(V b G,V a G)) K 2N 

= (X7 b F, A V a G) R « - (X7 a F, A T V 6 G), 



(7) 



where F,G £ C 1 (M) and where V& and V a denote the gradients with respect 
to b = (pi, . . . , 6jv) and a = (ai, . . . , ajsr), respectively. Therefore, equations JSJ 
can alternatively be written as b n = {b n , H}j, a n = {a n , H}j (1 < n < N). 

Since the matrix A, defined by iJBJ, has rank N — 1, the Poisson structure J 
is degenerate. It admits the two Casimir functions 2 



Ci := -— &„ and C 2 := I J| a„ J 

ra=l \n=l / 



(8) 



The gradients V&,aCj = (V&Cj, V a C,) (i = 1, 2) are given by 



1 



V*A = --(l,...,l), v a c 1= o, 

V b C 2 = 0, V a C 2 = ^ ( i . . . , — ] 
W Van ajv, 



(9) 
(10) 



They are linearly independent at each point (b, a) G M. 

Let M/3,c := {(&, a) e M : (d, C 2 ) = (/?, a)} denote the level set of (Ci, C 2 ) 
at (/3, a) e M x M>o. As Ci and C 2 are real analytic on M and the gradients 
Vh, a Ci and Vf, ja C 2 are linearly independent everywhere on M, the sets JAp, a 
are real analytic submanifolds of M of (real) codimension two. Furthermore the 
Poisson structure J, restricted to M/3 iCn becomes nondegenerate everywhere on 
and therefore induces a symplectic structure on M^q,. In this way, we 
obtain a symplectic foliation of M with M^ Q being the symplectic leaves. We 
denote by Hp a the restriction of the Hamiltonian H^oda to a . 

Besides Flaschka's coordinates of the Toda lattice we will also need to con- 
sider relative coordinates. Introduce . . . , vn) G given by Vi := qi+i — qi 
for 1 < i < N — 1 and ujv := ^ £i=i ft- Then u = Mq is a linear change of the 



2 A smooth function C : M — » R is a Casimir function for J if {C, -}j = 0. 



■5 



position coordinates where M is the N x iV-matrix 

/ -1 1 ... \ 

'•■ '■• : 

M - . 
...-11 

ViV- 1 N- 1 J 

The variables u = (u\, . . . , u^) conjugate to v — (vi, . . . , vn) are then given 



by u — (M 1 ) 1 p. The matrix (M T ) 1 can be computed to be 



(M 



T\-l 



1 

N 



/I 

2 



ViV... 



1\ 



/ 1 o 



11 ... 
1 ... 10 

\o o/ 



For any {(3, a) G K x M>o the variables (ffc, Ufc)i< / t<Ar_ 1 are canonically 
conjugate variables on M^ jQ . We want to express the Hamiltonian Hp_ a in 
terms of the coordinates (v,u) — {v k ,u k )i< k <N-i- Note that on M/3 iQ , u k = 
k(3 — Y^j=\Pj for 1 < fc < iV — 1 and um — N(3. Hence p\ = —Ui + (3, 
Pn = un-i + (3, and p k = {u k -\ — u k ) + (3 for 2 < k < N — 1, and thus 

1 V > 2 N [3 1 / n / \2 / \2 2 \ 

x2^Pj = + o ( U l + ("1 " U 2) +... + (MjV-2-«AT-l) +M Ar _ 1 ). 



Moreover, using that - gjv+i = qw - qi = Ylk=i(lk+i - 9fe) one gets 



N 



JV-1 



3=1 



Combining the two expressions displayed above yields 



Hp tCe : 



NP 2 1 2 



N-2 



/N-l 



+^ U?+E^-"w) 2 +^-i W E ^ +^ H . (ii) 



\k=i 



The coordinates (v k , u k )i<k<N-i (as well as u N , but not u^v) can easily be 
expressed in terms of (bj, dj)i<j< N} 



N 



w fe = 21og— , u k 

ak 



N 



E^ + E 6 ^ 0-<k<N-l) 



(12) 



3 = 1 



3 = 1 



Note that u k depends linearly on the 6-coordinates and is independent of the 
a-coordinates. On the other hand, v k depends only on the a-coordinates. The 
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partial derivatives of the Vk's at a = ctljv, where ljv = (1, • • • , 1) € l w , can be 
computed to be 

: 6 v + ( 13 ) 



daj a Na 



3 Birkhoff coordinates 

Let us recall from ^ and |S] our results concerning global Birkhoff coordinates 
on the phase space M. As a model space, we introduced the space CP := M. 2N ~ 2 x 
M x M >0 , foliated by V/3, a '■= R 2Ar ~ 2 x {/3} x {a} which are endowed with the 
standard symplectic structure. Denote by Jo the degenerate Poisson structure 
on CP having ?0, a as its symplectic leaves and the coordinates (3 and a as its 
Casimirs. 

In 0] we have proved that the Toda lattice admits global Birkhoff coordi- 
nates. 

Theorem 3.1. There exists a map 

$ : (M, J) -> (CP, Jo) 

(b,a) i-> (OEfc,yfc)i<fc<jv-i,Ci,C 2 ) 

™£/i £/ie following properties: 

• $ is a reaZ analytic diffeomorphism. 

• <1> is canonical, i.e. it preserves the Poisson brackets. In particular, the 
symplectic foliation of M &?/ M/3 jQ is trivial. 

• TTie coordinates {xk,yk)i<k<N-i,C\,C2 are global Birkhoff coordinates 
for the periodic Toda lattice, i.e. the transformed Toda Hamiltonian H* = 
iJo$ _1 is a function of the actions (Ik)i<k<N-i aniC\,Ci alone, where 
Ik = \{x 2 k + y 2 ) f or an V ^ < k < N — 1. (In the sequel we will denote H* 
by H as well.) 

As an immediate consequence of Theorem 13.11 one gets 

Corollary 3.2. For every (3 £ K and a > 0, 

$(M/j, ) = fp, a , 

and $\m/3 Q : 3Vtg,a — > CP^.q, is a symplectomorphism. In particular, the moment 
map 

is onio. 

We wish to analyze the Birkhoff map $ and its restriction to the leaves M,p, a 

near the equilibrium points (b, a) = (ftljsr, ctljsr), where ljy = (1, , 1) £ R w . 

Introduce for k £ Z with 1 < |fc| < N — 1 the abbreviation 
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To compute the Jacobian of <&~ 1 at the points (Oat-i, Ojv-x, j3, a) we first 
compute the Jacobian of <&. It turns out that the formulas are easier to express 
in complex notation. In |2] we have shown the following 

Lemma 3.3. For any (3 6 R, a > 0, and 1 < k < N — 1, the gradient 
^b,a( x k +iyk) at (b,a) = ((3l N ,al N ) is given by 



yb,a( x k + iyk)(j) = 
and 



1 



1 



3 2TTl(j-l)k/N (1 < j < N) 



(14) 



^2a~N A fe \ -2^~ 1 ) k ' N (N + 1 < j < 2N) 
V 6 ,„Ci = -4(ljv,0iv), V M C 2 = i(0 A r,l JV ). (15) 



According to the formulas (|14|) and (|15|) . the Jacobian (4^$ of $ at (6, a) 
(/31/v,alAr) is given by the (2N x 2iV)-matrix 



COS 



(2j-2)7rfc 
N 



V \/2aJV >fc 



sm 



(2j-2)7rfc 
N 



-N- 



N 



V o 







COS 



(2j-l)7Tfc 

N 



\ 



kj 



( 2_i 

V V2aJV <^fc 



sm 



(2j-l)xfc 
JV 



(16) 



where in each of the four (N — 1) x A^-submatrices the row and column indices 
k and j run from 1 to N — 1 and 1 to TV, respectively. In order to compute 
its inverse, we note that (|16f> can be written as the product A 2 • P ■ Ai of the 
diagonal matrices 

Ai := diag (ljv, -2 • ljv) 

and 



-Ldiag ( — 



1 



1 



1 



with the orthogonal 2N x 2N matrix 



w-i 






/ P (i) 


p(2)\ 


1 


p(3) 






In 







V o 


Ijv / 



where for 1 < i < 4, the submatrices (Pf^) are (Af — 1) x A^-matrices given 

(2j~2)7Tfc ^ ; p(2) - = ( cos (2j-l)^ 



byl^ := (cos x 

(2j-l)xfc \ 
JV /fcj 



JV 



fry 



, P( 3 ) := 



sin 



(2j-2)7Tfc> 



p(4) . = 



sm ■ 
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The inverse (d& i0 $) 1 = A± 1 P T A 2 1 at (6, a) = (j3 ■ In, a ■ In) can then 
easily computed to be 



/ 



^A.cos^M 



. V Afe cos — iv — 



J* 



V 



A fc sin 



(2j-2)7rfc 
N 



Ok 



V 2JV Afc Sm JV 



3* 



-10\ 







(17) 



again with 1 < j < JV and 1 < < AT — 1 , but with fc and j now being column 
and row indices, respectively. 

The row vectors in the 2iV x 2 TV-matrix 1)16(1 are the gradients Vb, a Xk, V^aZ/fc 
(1 < k < N — 1) and Vb, a Ci (i = 1,2). In the sequel we will consider the 
restrictions of the coordinates (xk,Vk)i<k<N-i to the symplectic leaves M.p >a . In 
order to obtain the gradients of these restricted functionals we need to take the 
orthogonal projections of the gradients Vb,a x k> Vb,aVki ^b,aCi onto the tangent 
space X^M^a of the leaf Mp, a at (b,a). Note that T h ^p, a {= R 2 ( JV ^ 1 )) is 
the orthogonal complement of span(Vt, !a Ci, Vf^aC^)- In view of the formulas 
(f^Tf) — (fTTTf) for the gradients of C\ and C2 one gets 



T b , a Mi3, a = {^,ri) G 



*2N 



2> = °.E£ 



0}. 



(18) 



3=1 



3=1 



For (b, a) = (/31jv, Qfljv); the conditions in the set in (|18f) simply are X^=i £7 = 



and X)i=i ^3 = 0- By formula (|14|l . one sees that for 1 < k < N — 1, 



3=1 J 



1 1 N 

(x k +iyk)= , t V e 27r 

%/2^V Xk ^ 



i(j-l)k/N _ 



3=1 



~i da 



— (x k + iyk) = - 



1 



V2aN A fc 



N 

e iKk/N e 2Mi-i)k/N 



0. 



Hence for (b,a) = (/31jv, aljv) the gradients V&, a ^fc, ^b,aUk are contained in 
T^aJVC^Q, for any 1 < < jV — 1, and there is no need to take projections. 



4 Linearized Birkhoff Coordinates 

For any given values of the Casimir functions, C\ = /3, C% = a, we wish to com- 
pute the first few coefficients of the Birkhoff normal form of the Toda Hamilto- 
nian near the elliptic fixed point (x, y) — (0, 0). We recall that these coefficients 
are essentially unique, i.e. do not depend on the choice of the Birkhoff coor- 
dinates. A possible way of proceeding is to substitute (b, a) = $ _1 (x, y, (3, a) 
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into the expression for the Toda Hamiltonian H = i Ylj=i &j + SjLi a ] an d 
then expand <& _1 (x, y, /?, a) at (x,y) = (0,0). However it seems difficult to 
explicitly compute the terms of this expansion, except the first ones which we 
have computed in the last section - see formula B17JI . Wc proceed differently 
and choose as the starting point of our computations the canonical coordinates 
Vk, Uk introduced in section rather than the non-canonical variables bj, aj. 
When expressed in these coordinates, the Toda Hamiltonian takes the form 
Hp, a = 2f- + H u + a 2 H v where by fTTJl. 

H u = J2(ui -u l+1 ) 2 + u 2 N ^ , (19) 

N-l 

H v = J2 e ~ V ' + e E ™^T (20) 
i=i 

Note that the Taylor expansion of Hp^ a at (v, u) = (Ojv-i, Ojv-i) is not in 
Birkhoff normal form up to order 2. In a first step we therefore want to choose 
a linear canonical transformation r)k)i<k<N-i l_ * ("^fe, Ufc)i<fe<Ar-i so that 
when expressed in the new variables (£,77) = r)k)i<k<N-i, the Toda Hamil- 
tonian is in Birkhoff normal form up to order 2. Consider the composition 

Up,* : R 2N - 2 (= Vp, a ) -► Mp, a -► R 2N - 2 

(xk,yk)l<k<N-l !-> {bj,dj)i<j<N !-> (vk,Uk)l<k<N-l 

of the inverse of the Birkhoff map ($|m^ Q ) 1 : K. 2Ar ~ 2 — > M.p >a with the coor- 
dinate transformation defined in (|12Jl . Then 0/3. Q is a canonical real analytic 
transformation as both (x,y) and (v,u) are canonical coordiantes for M/3 jCe . Its 
Jacobian 

R , a : R 2N ~ 2 -> M 2JV - 2 , g,»))w(t;, U ) = 4Aa|(x,y)=(o,o)(C 1 r)) (22) 

at (a;, y) = (0, 0) is a linear transformation with the desired properties. We will 
compute i?/3,ct as a composition of the Jacobian of (x, y, (3, a) 1— » (6, a) at (x, y) = 
(0,0) (with /?, a fixed) and the one of (6, a) 1— > (v,u) at (6, a) = ((31n, aljv). 
It is convenient to use complex notation for £,k, Vk (I < k < N — 1), 

Cfe == ~^(^ k ~ := ^7f ^ fc + i?7fc )' 

where the sign in the definition of is chosen so that dC,k A <i£_fc = *c££jfc A drjk- 
The vector £ = (Cfc)i<|fc|<iv-i is an element in the space 

Z := {C = (Cfc)i<|fe|<JV-i G C 2Ar - 2 : C-fc = Cfe" Vl<fc<7V-l}. (23) 
The components of C satisfy the identity 

e W Cfc + e -*tf*/" C _ fc = ^ (cos £fe + sin % ) . (24) 
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As for any 1 < k < N — 1, Uk is a linear function of the frj's by (|12fl one has 

N k 



k 

3 = 1 3 = 1 



where 6j(C) (1 < j ' < N) can be computed using formula (|24[) and the Jacobian 
of ($1^, J^ 1 at (x,y) = (0,0) obtained in JT7J) to get, for 1 <j <N, 

JV-i 

6,(C) = E A ^ (e^'-^Cfc + e-^-^^C-fc) 

k=i 

l<|fc|<iV-l 

where we used that A^ = for any 1 < k < N — 1. Note that J2j=i ^j'(C) = 
as Ef=i e^-D*/" = 0. To compute H u = \{u\ + E^ 2 ~^+i) 2 + 
we use 

ui(C) =6i(C) - ( 25 ) 

1<|A|<JV-1 

u,(C)-« I+ i(0=-6l+i(0 = -V^V : E ^e^ N Ck, (26) 

l<|fc|<AT-l 

-UN-xiQ = - E & i(0 = MO = £ A fc e 2 ' r< ( w - 1 ) fc / Jv a-(27) 

J=l l<|fc|<JV-l 

with 1 < Z < JV- 2 in fZ®. Thus 

^(o = ^e( E A fc e 2 ^^a 

i=0 \l<|fc|<iV-l 
= ^ E A fe A fe ,fx;e 2 ^)/^ aCfc ,. 

l<|fc|,|fc'|<iV-l \ i=o ) 

Using that Ya=o e 2nilk / N = N5 k0 for any < |fc| < N - 1, one gets an 
expression which is quadratic in the ^-variables, 

N-l 

ffn(C)=«E A 'CfcC^ fc . (28) 

k=l 

Now let us turn to given by formula J20J). As (6, a) = (/31jv, aljv) we 
have by ljP3)l. 

2 ... 2 A ... 
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where a,j(() (1 < j < N) can be computed using formula 1(24(1 and the Jacobian 
<fc,j/,/3,a$ _1 at (x , y) — (0,0) obtained in l(T7)l to get 



JV-l 

k (j*(2j-l)k/N ^ + e - l7r ( 2j -l)fc/JV c _ fc 

fc=l 



fe=i 



l<|fe|<iV-l 

Again we have that Y2j=i a j(0 — as for any 1 < |fc| < N — 1, 

JV / N 

y^ e m{2j-l)k/N = e -iirk/N y^^ijk/N I = q 

i=l \J=1 

Hence, for 1 < i < JV-1, 

«; = ^= V A^e^e-^/^a- (29) 

v l<|fe|<iV-l 

Define «o by the expression on the right hand side of 1(29(1 evaluated at I = 0. 
Note that 

JV-l /N-l \ 

E vi = -i E A fe C fe e— fe / w E e2 ™ ifc/JV = °- ( 3 °) 



1=0 v Ltx ' l<|fc|<JV-l \ i=l 

Hence Ez^i 1 u ' ~ ~ w o and therefore 

ff« = E e ~ 

Z=0 

Now we expand H v at w = up to order 4, 

N-l „ AT-1 „ jV-1 , N-l 



N-l 



Hv=N -Y,m + \Y.$- t\ E^ 3 + jy Etf+W)- ( 31 ) 



2^ 1 3! ^ 1 4! 

=0 2=0 2=0 /=0 



By JSIll, EzLo «i = 0. Substituting the expressions 1)29(1 for u; in the quadratic 
term in the expansion 1)31(1 we get 



N-l , /N-l 



2=0 i<|fc|,|fc'l<^v-i \ J =° 

2 

= - E x l(kC-k, 

fe=l 



where we again used that Afe = A_fc and that ^"q 1 e 2^Uk/N _ jv,5 fc0 f or an y 
< | A; | < AT — 1. The terms of third and fourth order in H v are treated similarly. 
Combining the above computations leads to 
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Proposition 4.1. Let (3 € K and a > 0. TTien 

ff/3,a o Rp, a {0 = Go + aG 2 + V^G 3 + G 4 + 0(C 5 ) (32) 

where Rp t a is the linear canonical transformation introduced in \2ty and Gi 
(t) < i < A) are given by 

G :=^ + Na 2 , (33) 

AT-l 

G 2 := 2 ^ AfoC-fc, (34) 

k=l 

G 3 :=--4= E (-l) (fe+fc ' +fc " )/Ar A fe A fe ,A fe ,C fc a< fc », (35) 

1 < I fc I , I fc "|<JV-1 
Hfe' + fe" = 0"iodJV 

G4:= 2i7V £ (-l) (fe ^"^"' )/Ar A fe A fe ,A^A^a.a<^C^. (36) 

l<|fc| ) |fe'| ) |fe // | ) |fe /// |<JV — 1 
k + k' + k" + k'" = mod N 

Note that Hp, a o R/3. a (C) depends on (3 only through the constant term 
and that it is in Birkhoff normal form up to order 2. 

To finish this section let us express the Birkhoff coordinates (x, y) in terms 
of the coordinates (£, rf) near the origin. The two coordinate systems are related 

by 

(x, y) = (ilp, a ) 1 do,ofy8,a(£,»7) 
where we used that Rp t a = do,o^/3,a- Hence 

(x k ,yk) = ^k,Vk) + 0(\(^r,)\ 2 ) Vl<k<N-l (37) 

and 

/ = d 

2 2 

Denote by C^ a the Poisson algebra of Hp. a o i?/3. Q , i.e. the space of germs 



Ik = 4 + rt = ti + 4 +omv)l 3 ) vi< k <N-i. 



of real analytic functions F(£, n) = J2\ 7 \+\s\>2 fisCrf sucn that {F, Ik} = for 
any 1 < A; < N — 1. In view of l|37|) we say that C£ a is non-resonant. The 
following result is then well known (see e.g. |7j, Appendix G). 

Corollary 4.2. For any f3 S R, a > 0, and to > 3, t/iere exists a (germ 
of a) real analytic canonical transformation of the form "Id + higher order 
terms" given by X F \t=i> where X F is the flow of the Hamiltonian vector field 
associated to the Hamiltonian 



7 Vi^ 

/70S 

3<l7l + l«5|<m 



such that F o X F \ t —i is in Birkhoff normal form up to order to for any F in 
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5 Proof of Theorem 11.11 

We now transform the Hamiltonian -ff/3 lQ of?^ ja (£) into its Birkhoff normal form 
up to order 4 by a standard procedure - see e.g. section 14 in 7]. The phase 
space Z, defined in l|23|) . is endowed with the Poisson bracket 

{F, G\ = i 2^ a k 



d( k dC-k ' 

l<|fc|<JV-l 

where = sgn (k) is the sign of k. The Hamiltonian vector field Xp associated 
to the Hamiltonian F is then given by Xp — i Y^i<\k\<N-i CTfc ;^^af^- With a 
first canonical transformation we want to eliminate the third order term G 3 in 
the expansion (|32|1 of if/3, a, o i? ( g Q ,(C). We construct such a canonical transfor- 
mation on the phase space Z as the time-l-map ^1 := X f F \ t= i of the flow X F of 
a real analytic Hamiltonian F := a~ 1 ^ 2 F 3 which is a homogeneous polynomial 
m (k (1 < \k\ < iV — 1) of degree 3 and solves the following homological equation 

{aG 2 ,a-?F 3 } + a^G 3 = 0. (38) 

To simplify notation we momentarily write H instead of Hp. a o Rp^ a . Assuming 
for the moment that Ij38(l can be solved and that X F is defined for < t < 1 in 
some neighbourhood of the origin in Z we can use Taylor's formula to expand 
H o X F around t = 0, 

HoX F = HoX%+ [ ^-{H o X F )ds = H + [ {H,F}oX s F ds 
Jo ds J 

= H + t{H,F}+ [ (t-s){{H,F},F}oX s F ds. (39) 
Jo 

When evaluating this expression at t = 1, one gets 

H o = G + aG 2 + {«G 2 , F} + [ (1 - t){{aG 2 , F}, F} o X%dt 

Jo 

+v^G 3 + / {V^G 3 ,F}oX F dt + G 4 + 0{( 5 ). 
Jo 

Using that {aG 2 ,F} + \/aG 3 = 0, the latter expression simplifies and we get 

HoV 1 =G + aG 2 + [ t{V^G 3 ,F}oX F dt + Gi + 0(C 5 ). 
Jo 

Integrating by parts once more and taking into account that F = a~ x / 2 F 3 is 
homogeneous of degree 3 one obtains, in view of (|39|l . 

Hp, a o R^ a o = Go + aG 2 + G 4 + ^{G 3 , F 3 } + 0(( 5 ). (40) 
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Note that {G3, F3} is homogeneous of order 4. Hence our first step is achieved. 
It remains to solve (|38ll . 

By Corollary 14.21 with m = 3 there exists a polynomial homogeneous of 
degree 3, 

l<|fe|,|fe'|,|fe"|<AT-l 

(with W kk , k ,, invariant under permutations of A;, A:', fe ) so that the time-l-map 
-X"v^|t=i of the flow X^y corresponding to the Hamiltonian W = a~ x / 2 W 3 brings 
any Hamiltonian in Co a into Birkhoff normal form up to order 3. In particular, 
the identity {G2, W 3 } + G3 = is satisfied. Note that 

{G 2 ,W 3 } = t Yl ^kX 2 k C- k m '' 



"dC-k 

l<|fc|<iV-l s 

-i Y ( Sfc + Sk ' + s k")wjf k } k „(k(k>(k" 

l<|fe|,|fe'|,|fc"|<JV-l 



as Sk = 2 sin ^ — 2<7fcAjj, and it follows that 



(s k +sv+s k «)iW$ k ,,=GW k ,, (41) 



where G k 3 k , k „ are the coefficients of G3 defined by lp?5|) . 

G (3) = f-l/(6\/iV)(-l)( fe + fe '+ fe ")/ A 'AfcAfc'A fe - forfc+fc' + fc" = 0mod7V 
kk'k" 1 q otherwise. 

As G { kk \ k „ ^ iff k + k' + k" = mod N it follows 3 from gTJl that for any triple 
(fc, fc', A") with k + k' + k" = mod TV, 

Sfc + Sfe' + Sfe» 7^ 

and 

)W (3) _ U kk'k" 

lYV kk'k" 



Sk + Sk' + Sfc" 

p (3) 

l<|fc|,|fc'|,|fe"|<JV-l r kk'k" 



Now define F 3 = Ei<ife|.|fe'|,|fc"|<jv-i F kk' ' k" CfeCfe' Cfc" and F = a VSjPg by 



F (3) f = f W^,, for k + k' + k" = mod N ^ 
kk ' k " \ otherwise. 

Then clearly {G2,F 3 } + G3 = and hence by the considerations above "Ji := 
Xp\ t= i has the property that Hp iCt o Rp^ a o satisfies identity J3QJ. 

Now let us investigate the 4th order term G 4 + \{G 3 ,F 3 } in (|3U|) . We 
decompose this sum into its contribution to the Birkhoff normal form and the 



3 See II 1 1 (cf. also 0) for a direct proof of this nonresonance condition. 
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rest, to be transformed away in a moment. Let us first compute {G3, F$} in a 
more explicit form. By (|35|l and l|42l) . 

dG 3 8F 3 



{G 3 ,F 3 } = * Y, 



Ok 



l<|fc|<AT-l 



36N 



- E 

l<|fc|<iV-l 
/ 



<7fc 



3 J2 (-l) r AfcA/A m OCr, 



i<|i|,| m |<jv-i, 

Z + m=-fc + r-JV 



3 ^ (-i)-' _ XkXi ; Xmi 



l<|i'|,|m'|<JV-l, 



S_fe + S// + S T , 

r + r ' SfeA;A m A/'A r 



K|fc|<JV-l l<|l|,|m|,|l'|,|m'|<N-l 
i+m-rW=-fc 
l' +m' -r' N = k 



S-k + Si' + s, 



~ ClCmCl' Cm' > 



where for the latter equality we used again that 2a k )? k = s k . Setting 

l + m + l' + m' 



£ll7ll'' 



N 



and using that s-k = — Sfe one then gets 



1< fe <iV-l I+m=-k™od» 



A; A m A;' A TO ' 



{G 3 ,^ 3 }-^ E 
i<|fe|</ 

iV-l 

iN At . A .. v ' -l + {s v +s m ,) s k 



-1 + (s v + s m i)/sk 
A;A m A;/ A m / 



CiCmO'Cn 



OCmO'C" 



fc — 1 l + m=-kmodN 
l'+m' = k mod N 



N-l 



^E E M) e 



A;A m A;' A 



87V 



fe=l l + m = kmod N 

l'+m' = — k mod N 



-1 - (sj/ + s m ,)/s k 



ClCmCl'Cn 



iV-l 



87v E E 

fc=l l+m=-tmoJ» 
i' + m' = fc mod 



1 



1 



-1 + (s// + S m ')/ s * -1 - (s; + s m )/s fc 



• ( — l) £ ' ml ' m ' A;A m A/'A m 'CiCmO'Cm' 

where . . . stand for 1 < \l\, \m\, \l'\, \m'\ < N-l. Note that for k = V + m' + r'N 
with 1 < k < N — 1 and r' € Z we have 



Sfc = |s;'+m'|- 



Hence 



^{G3,F 3 } — — — ^2 Clml'm>{— l) Slrr >>' m ' XlX m Xl'X m 'ClCmCl'Cm' (43) 



24iV 

Z + m + i' + m'iEO mod N 
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where 4 



QmZ' 




l<|l|,|m|,|I'|,|m'|<JV-l 

I + to ^ mod N 

l + m + l' + to' = mod JV 



otherwise 

(44) 

Combined with formula ((36(1 for G4, the sum G4 + ^{Gs,^} equals 
24]v 51 (-l) £ "' fc "' t "'(l + cfefe'fe"fc»0 ■ AfcAfc'Afc»A fe /»aCfc'Cfc''Cfc'»- (45) 

1< |fc|, I ,| A;" , I A;'" |<7V— 1 
k + k' + k" + k'" = Q mod N 

We now decompose ((45(1 into its contribution itn(G4 + ^{G 3 , -F3}) to the Birkhoff 
normal form of Hp a o Rp^ a and the rest. 

Lemma 5.1. The normal form part of G4 + ^{G^jF^} is given by 

/ 1 \ 1 N ~ 1 

n N G4 + -{G 3) F 3 } = — ^ |0 1 4 - (46) 



4iV 

i=i 

Proof. The indices ft, ft', ft", ft'" of the terms in 145(1 contributing to the normal 
form satisfy the condition 

{ft, ft', k",k"'} = {l,-l, to, -m} (47) 

with l<£<m.<AT — 1. In the case I — to, {I, —I, I, —1} is viewed as a set-like 
object whose two elements / and — I each have multiplicity two. 

We investigate ttn(G4) and 7rjv(J; {G3, F 3 }) separately. Let us start with G4. 
We distinguish the cases I = m and I ^= to in (147(1 . For I — to, there are ( 4 ) = 6 
distinct combinations of indices (k, k' , k" , ft"') satisfying 1(47(1 . whereas for I ^ to, 
there are 4! = 24 (automatically distinct) permutations of (I, to, — I, —to). Hence 
we have 

N-l 



MG 4 ) = ^ 6]TAf|0| 4 + 24 A ? A ™K/| 2 |C 

y 1=1 l<Krn<N-l 

= 4^ E A ^i 4 +4 £ «Jc*i 2 ic m 

y i=l l<;<m<N-l 



(48) 



Now let us compute kn(\{Gzi F 3 }). We have to single out the matches of 
<(47|) for which in addition Ckk'k"k"' 7^ 0, i.e. 

k + k' ^ OmodTV, ft + ft' + ft" + ft'" = OmodiV 



4 To keep the formulas as simple as possible we have not symmetrized the coefficients 

^Iml 1 m' • 
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To fulfill H47fl . there are therefore the two possibilities 

k + k" = k + k'" = 

k' + k'" = OT k' + k" = ' ( > 

In both cases, we have Sk" + s/c'" = — (sfc + Sk>), and therefore l|44|) reduces to 

-3\sk+k>\ /f.p.x 

Ckk'k"k"' = i j— ; ; • \0V) 

\Sk+k<\ + Sk + Sk' 

Note that l|5l)|l remains valid for k + k' = N, since in this case Sk + k' = and 
Sk + s^ > as k and kl must satisfy 1 < k, k' < N — 1, but not for k + k! = 0, 
since in this case \sk+k' \ + Sk + Sfe' = 0. 

We first compute the diagonal part of itn(^{G3, ^3})- In this case, the two 
possibilities in (|49|l coincide and the solutions are 

^' k "^ = {(-i-l,~h~l)- (51) 

The sum of the coefficients Ckk'k"k"' for the two cases listed in (|51ll is 

■ . f 1 1 \ —6s?, „ 7 In 

cu-i-i + c-i _ u ,i = -3 s 2 z — h : — — = -= f-j = 6cot T7- 

\l s 2;|+2si \s2i\-2si) s^-Asf N 

We now turn to the off-diagonal part of ttn (^{G^, F3}) . The matches of 
(32} , 63) for given {I, to} C {1, . . . , JV — 1} with Z < to, (fc, fe') = (±Z, ±m) and 
(&",&'") = (±Z,±to) are 

(Z, m, — Z, —to) 
(Z, —to, — Z, to) 



^/c, k ^ k ^ k ^ — ^ 



(— Z, m, Z, —to) ' 



^ (-Z, -to, Z, to) 

The remaining matches are obtained from (|52|1 by permuting the first and second 
or the third and fourth columns on the right hand side of l|52(l , bringing the total 
of all matches to 16 = 4-4. Note that by formula QMljl. these permutations leave 
the value of the coefficients Ckk'k"k"' invariant. Taking the sum of the coefficients 
Ckk'k"k"' for all the matches, we obtain from (|50|l 

4(c/ ; m, — I, — m Cl,—17i,—l,m ~t~ C—l,m,l, — m C—l,—m,l,7n) 

_ _ 12 . |sj+m| 1 | S Z— ml 



i|+Si+S m |Sj_ m |+Sj 
-m| |s;+ m | 



\Sl-rn\ ~ S t + S m \Sl +m \—Sl-S 

= -24 



2 2 
l—m 1 Z + m 



= -24(2a?_ m a? +m - s^Jfj + f™g ~ s f+nM ~ *m) 2 ) 

S ?-m S r+m + {si-S m ) 2 {si+S m ) 2 -S^{sl + S m ) 2 -sl +m {si-S m )' 2 

= -24, 
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*l—m"l+m 



(si — s m ) 2 (si + s m ) 2 . Collecting terms, we thus have 



kn[\{G 3 ,F 3 }\ 



V 1=1 l<l<m<N—l 
y /=1 l<;<m<Ar-l y 



(53) 



Adding up (gSJl and (JHSJ, we obtain |gSJ|. □ 

In a next step we want to remove (1 — ttn) (G4 + \{G 3 , -F3}) from the Hamil- 
tonian by a second coordinate transformation \&2- In view of Corollarv l4. 21 with 
m = 4, there exists a canonical transformation "J 2 of the form "Id + higher 
order terms" so that for any Hamiltonian F in the Poisson algebra ^\C^ a , 
Fo$ 2 is in Birkhoff normal form up to order 4. We have proved the following 

Proposition 5.2. Let [3 £ K and a > be given. The real analytic symplectic 
coordinate transformation £ = S(z) = fio ^2 (2) defined in a neighborhood of 
the origin in Z, transforms the Hamiltonian Hp i0 o Rp a into its Birkhoff normal 
form up to order 4. More precisely, 



H , a o R p , a o H = G + aG 2 + n N (g a + i{G 3 , F 3 }^j + 0(z 5 ), 



(54) 



with Gq, Gi and ttn (G 4 + ±{G 3 ,F 3 }) defined by and g3J). 

Proof of Theorem\Tl\ The map ty? )Q : 7 p>a -> K 2Ar ~ 2 introduced in P is a 
canonical real analytic diffeomorphism so that Hp a o fin „ — (Ha a o Rp a ) o 
(JJ^ a o £lp ia ) is in Birkhoff normal form. By the definition (|22|) . the map Rp a ° 
£lp, a is canonical and, near the origin, of the form "Id + higher order terms". 
The canonical transformation 5 of Proposition 15.21 is also of the form "Id + 
higher order terms" and brings Hp^ a oRp a locally around the origin into Birkhoff 
normal form up to order 4. Hence the transformations Rp \ ° Qp, a and 3 differ, 
near the origin, by a transformation of the form "Id + higher order terms" . 
By the uniqueness of the Birkhoff normal form, the expansions of the Toda 
Hamiltonian near the origin, when expressed in these two coordinate systems, 
coincide up to order 4. Hence Proposition l5 . 21 provides us with the Taylor series 
expansion of < Kp, a (I) — — h H a (I) in terms of the actions 

t — it \ t — x k + y"i 

1 — \J-k)l<k<N-li 1 k — ^ 

up to order 2. In view of 133|) . I|34|) . and l|46|) one has 

N-i N-i 

k=l k=l 

This proves Theorem ll.il □ 



REFERENCES 



19 



References 

[1] G. P. Berman & F. M. Izrailev, The Fermi-Pasta-Ulam problem: 50 
years of progress. arXiv:nlin. CD/04 11062 v3 

[2] H. FLASCHKA, The Toda lattice. I. Existence of integrals. Phys. Rev., Sect. 
B 9 (1974), 1924-1925. 

[3] M. Henon, Integrals of the Toda lattice. Phys. Rev., Sect. B 9 (1974), 
1921. 

[4] A. Henrici & T. Kappeler, Global action-angle variables for the periodic 
Toda lattice. Preprint, 2006. 

[5] A. Henrici & T. Kappeler, Global Birkhoff coordinates for the periodic 
Toda lattice. Preprint, 2006. 

[6] A. Henrici & T. Kappeler, Birkhoff normal form for periodic FPU 
chains. Preprint, 2006. 

[7] T. Kappeler & J. Poschel, KdV & KAM. Ergebnisse der Mathematik, 
3. Folgc, 45, Springer, Berlin, 2003. 

[8] S. V. MANAKOV, Complete integrability and stochastization of discrete dy- 
namical systems. Zh. Exp. Teor. Fiz. 67 (1974), 543-555 [Russian]. English 
translation: Sov. Phys. JETP 40 (1975), 269-274. 

[9] J. Poschel, Integrability of Hamiltonian systems on Cantor sets. Comm. 
Pure Appl. Math. 35 (1982), 653-695. 

[10] J. Poschel, On Nekhoroshev's Estimate at an Elliptic Equilibrium. Int. 
Math. Res. Not. 4 (1999), 203-215. 

[11] B. Rink, Symmetry and resonance in periodic FPU chains. Comm. Math. 
Phys. 218 (2001), 665-685. 

[12] M. Toda, Theory of Nonlinear Lattices, 2nd enl. ed., Springer Series in 
Solid-State Sciences 20, Springer, Berlin, 1989. 



